## Example of a primitive polynomial

from PyM import *


K = Zn(3)

# Ring K[T]
[_,T] = polynomial_ring(K,'T')

f = T**4 + T - 1

Q = is_irreducible(f,K)

show(Q)

[F,x] = extension(K,f,'x')

show(order(x))

show(primitive_root(F))